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Abstract 

With the help of the Diffusion Entropy technique we show the non-Poisson statistics 
of the distances between consecutive Omori's swarms of earthquakes. We give an 
analytical proof of the numerical results of an earlier paper [Mega et al., Phys. Rev. 
Lett. 90 (2003) 188501]. 
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The model usually adopted to describe the time distribution of earthquakes is 
the Generalized Poisson (GP) model [1,2,3,4,5]. The GP model assumes that 
the earthquakes are grouped into temporal clusters of events and these clusters 
are uncorrelated, and, therefore, completely unpredictable: The clusters are 
supposed to be distributed at random in time and therefore the time intervals 
between one cluster and the next one follow a Poisson distribution. The intra- 
cluster earthquakes are in fact correlated as it is expressed by the Omori's law 
[6,7], an empirical law stating that the main shock, i.e. the highest magnitude 
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earthquake of the cluster, occurring at time to is followed by a swarm of 
triggered earthquakes (after shocks) whose number (or frequency) n{t) decays 
in time as a power law, n{t) oc (t — to)~'', with the exponent p being very close 
to 1. If we denote with the symbol r the time intervals between one earthquake 
and the next, then right after a main shock, a short value of r is followed with 
a large probability by another short value. For the same reason, far from a 
main shock and prior to the next one, a long value of r is followed by another 
long value of r. This implies that the correlation function ((rj — {T)){Tj — (r))) 
is not zero for i j and that it survives for all pairs of seismic events in 
between two consecutive unpredictable shocks. Omori's law also implies [8] 
that the distribution of r, is a power law ip{r) oc r^^. This espression is valid 
in the time regime inside a swarm, and it is then truncated by a sharp cutoff 
caused by the arrival of the next swarm. 

The catalog we have studied covers the period 1976-2002 in the region of 
Southern California spanning 20° N -45° N latitude and 100° W 125° W lon- 
gitude [9]. This region is crossed by the most seismogenetic part of the San 
Andrea fault, which accommodates by displacement the primarily strike-slip 
motion between the North America and the Pacific plates, producing veloc- 
ities up to 47 mm/yr [10]. The total number of recorded earthquakes in the 
catalog is 383687. 

Herein we disprove the GP model, providing evidence for the non-Poisson 
statistics of inter-cluster times, by applying to the mentioned catalog the Dif- 
fusion Entropy (DE) technique [11]. Here we discuss with analytical arguments 

some issues that in an earlier paper [12] we have examined by means of a nu- 
merical treatment. As in the GP model, we assume that each cluster starts 
with an unpredictable triggering event (it may or may not be the main-shock) . 
The distance between the i-th and the i + 1-th cluster is therefore the time 
distance between such events, which we indicate as r]"*^, obeying the non- 
correlation property, (rl^'rj™^) = (r'™!)^ if i j. The superscript m stands 
for main-shock, but we actually need not to make the assumption that the trig- 
gering event is a large earthquake. As we shall see, the DE measures statistical 
properties of events with no need of identifying them. 

Let us recall the definition of the DE functional S{t) as the Shannon entropy of 
p{y, t), the probability distribution to observe a fixed number of seismic events 
y in a given time interval [11]. This observation is equivalent to observing the 
spreading of a number of walkers making one step forward at each time where 
an event is met. Hence the term "Diffusion Entropy". Different trajectories 
are chosen with the usual method of observing different time windows in the 
sequence, which is herein assumed to be stationary. Let us assume that p{y, t) 
follows the scaling law with the form 
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where 5 is a positive exponent and F{x) is a positive and integrable function 
of X. As a consequence of this assumption, after a straightforward algebra, we 
find that 

S{t) + 5 \n{t). (2) 



This means that the entropy of the diffusion process is a hnear function of ln(t) 
and a measure of the slope is equivalent to the determination of the scaling 
parameter S. 

We now show that, in the earthquake series under investigation, the principal 
source of entropy increase is given by the occurrence of the cluster-initiating 
seismic events. Let us indicate with 0(t["*1) the probability density function 
(pdf) of times between chistcrs and with h{x) the pdf of the number of earth- 
quakes in a cluster. The function h{x), usually referred to as the Pareto's law 
of earthquakes, is known to decay as h{x) ~ l/x°'^^, where a is a positive 
number. In the literature an exponent for clusters size distribution a ranging 
from a — 1.25 [13] to 3 [4] is reported. On the other hand, we assume for 
0(t["*1), which is exponential in the GP model, a form ^(rt™!) ~ (l/r'"*!)'*. We 
also assume that 2 < /i < 3. As it will become clear in the next paragraphs, 
beyond the upper limit it is impossible for the DE to distinguish between an 
inverse-power law and an exponential; for < 2, on the other hand, the signal 
cannot be stationary [11]. The connection between x and the early defined 
variable y, over which the DE is calculated, is 

zit) 

y{t) = Y.^r, (3) 



where the sum is carried over different clusters, Xi is the number of shocks in 
the i-ih cluster and z{t) is the number of clusters in the same time window of 
length t considered for y. 

Let us call 4>{s) the Laplace transform of (/)(r["*l) and h{k) the Fourier trans- 
form of h{x). Suppose now that the time duration of a cluster is negligible 
with respect to the mean time distance (r'™"!). In this case we can use directly 
a Continuous Time Random Walk (CTRW) formalism to calculate the prob- 
ability piy, t) that a random walker, moving of a quantity x at the time at 
which there is a cluster of size x (and resting otherwise), is at position y after 
a time t. For the Fourier- Laplace transform of t) we will have, from the 
theory of CTRW [14,15], 

p(k,s)^ — . (4) 
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To obtain the asymptotical behavior we write, if 2 < // < 3, 

${s) ~ 1 - (r[™l)s + cs^, (5) 



where j = /i — 1. Eq.(4) becomes 



p{k,s) 




(6) 



1 - h{k) + s(tH)^(A;) - ch{k)s'r 



Considering that 



h{k) ~ l + ik{x) + bk'', 



(7) 



if 1 < Q! < 2, we see that Eq. (6) leads to a ballistic scaling in the laboratory 
reference frame. The DE is insensitive to drifts, so we assume that Eq. (1) is 
fulfilled in a "detrended" moving reference frame, namely, where the position 
of the walker y fulfills the condition {y{t)) = at each time t. To obtain p{k, s) 
in such a moving reference frame we have to perform the following substitution 



After that, we apply the diffusive limit k ^ s obtaining 



thus proving that the most anomalous, namely, the smallest between the two 
exponents, either a + 1 or /i, determines the asymptotic scaling according to 
the prescription 5 — 1/a or S — l/{fj, — 1), respectively. 

In Fig. 1 we report the results of the DE method. In full circles we plot the 
entropy S{t) as a function of time when all the seismic events of the catalog are 
considered (independently of their magnitude M). A fit in the linear region 
gives a value of the scaling parameter 5 = 0.94. We next consider (open 
symbols in Fig. 1) only the earthquakes with magnitude larger than a fixed 
value M = 2,3,4. We see that, regardless of the value of the threshold M 
adopted, the function S{t) is characterized by the same long-time behavior 
with the same slope. This indicates that we are observing a property of the 
time location of the main earthquakes. Moreover, if the value of S were due 
to h{x) we should observe at most 6 = 1/a = 1/1.25 = 0.80. This leads us 
to conclude that the asymptotic form (and scaling) of p{y, t) is determined 
by the probability p{z, t) of finding z unpredictable events in a time window 
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Fig. 1. The Shannon entropy S{t) of the diffusion process as a function of the 
logarithm of time, expressed in minutes. From top to bottom, the curves refer to all 
events (full circles) and to events with threshold M = 2,3,4 (open symbols). The 
straight lines have the slope S = 0.94. 

of duration t. This is in turn determined by a 0(rt™]) decaying as an inverse 
power law with exponent = 2.06. Using the theory of Ref. [11] we also 
determine the form of p{z,t), which is well approximated by an asymmetric 
Levy distribution with Levy index // — 1. 

In conclusion, in this paper we have studied the statistical properties of earth- 
quakes time distribution. Inter-cluster distances obey an inverse power law 
prescription (/)(rH) oc (rH)-M with fi = 2.06 ± 0.01 thus ruhng out the GP 
model. The method proposed is based on the fact that the asymptotic proper- 
ties of diffusion process generated by the seismic events arc scarsely sensitive 
to the memory stemming from the Omori's law. They are, on the contrary, sen- 
sitive to the anomalous statistics generated by the non-Poisson nature of the 
time distance between two consecutive large earthquakes. This non-Poisson 
behavior reflects, in our opinion, the cooperative behavior of the geological 
processes triggering the main shock, and consequently, to some extent, some 
sort of predictability. The emergenge of this possibility will be investigated in 
the future. 

As a final remark, we recall that throughout our analysis we made the rea- 
sonable hypothesis that the series is stationary. On a formal ground, a scaling 
value 5 ~ 1 is compatible with a nonstationary process. These nonstationary 
contributions may have different orgins. The most trivial is the lack of statis- 
tics: we may have for instance, a small number of Omori's swarms, but very 
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extended in time. For this reason we repeated the analysis using only portions 
of the series, obtaining the same results. Another, more interesting, source of 
non-stationarity could be a strong correlation among the series of the trigger- 
ing processes, namely the condition (T|™''rj™'') 7^ (rt™!)^. It is interesting to 
notice that also this condition might reflect a form of predictability to assess 
through a properly tailored form of statistical analysis of time series. We leave 
this ClS Sb subject of further investigation. 
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